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Abstract 

An interacting scalar field with largish coupling to curvature can support a distinctive inflationary 
universe scenario. Previously this has been discussed for the Standard Model Higgs field, treated 
classically or in a leading log approximation. Here we investigate the quantum theory using 
renormalization group methods. In this model the running of both the effective Planck mass and 
the couplings is important. The cosmological predictions are consistent with existing WMAP5 
data, with 0.967 ns 0.98 (for Ne = 60) and negligible gravity waves. We find a relationship 
between the spectral index and the Higgs mass that is sharply varying for mh ~ 120 — 135 GeV 
(depending on the top mass); in the future, that relationship could be tested against data from 
PLANCK and LHC. We also comment briefiy on how similar dynamics might arise in more 
general settings, and discuss our assumptions from the effective field theory point of view. 
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1. Introduction 

The hypothesis that there was a period in the early history of the universe during which a local 
Lorentz invariant energy density - i.e., an effective cosmological term - dominated the equation of 
state, causing exponential expansion, explains several otherwise puzzling features of the present 
universe (flatness, isotropy, homogeneity) [T] |31 U]. It also suggests a mechanism whereby 
primordial density fluctuations arise through intrinsic fluctuations of quantum fields, leading 
to qualitative and semi-quantitative predictions that are consistent with recent observations. 
However the physics behind inflation remains mysterious. What, specifically, is the source of the 
energy density? Ideas ranging from fields associated with supersymmetry, string moduli, ghosts, 
branes, and others abound [HI El [3 HI IHl HH]- One (or more) of them might be correct, but all 
are highly speculative, and none is obviously compelling. 

Alternatively, we can look for inflationary dynamics based on degrees of freedom already 
present in the Standard Model. We can also attempt to maintain the guiding philosophy of the 
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Standard Model, including gravity, to allow only local interactions which are gauge invariant 
and have mass dimension < 4. Within this very restrictive framework, there remains the pos- 
sibility to include the non-minimal gravitational coupling ^H^ HTZ. Here H is the Higgs field, 
TZ is the Ricci scalar, and ^ is a dimensionless coupling constant, whose value is unknown and 
largely unconstrained by experiment]^ Indeed renormalization of the divergences arising in a 
self-interacting scalar theory in curved spacetime requires a term of this form [TT] . The Higgs 
sector is then described, classically, by the Lagrangian 

Ch = -\dH\'^ + ^?H^H - \{H^Hf + ^H^HTZ (1) 

where A is the Higgs self coupling and fi is the Higgs mass parameter. 

It has been known for some time that such minimal classical Lagrangians can support inflation 
driven by an interesting interplay between the quartic term and the non-minimal coupling term 
[HI [inj [in [T5]. For ease of reference, we will call this general set-up "running inflation"; the 
name seems appropriate, since evolution of the effective Planck mass and the effective scalar 
mass is central to the dynamics]^ This quasi-renormalizable set-up allows use of renormalization 
group methods, as will be illustrated here. By quasi-renormalizable, we mean that the theory is 
renormalizable when gravity is treated classically; in particular, we ignore quantum corrections 
from graviton exchange (see Appendix |b]) . In the investigation of (non-gravitational) quantum 
effects, it is appropriate to focus specifically on the Standard Model, for two reasons. First, 
because (as we will see) it illustrates important qualitative issues in a very concrete, familiar 
setting. Second, because it - or something close to it - might actually contain the degrees of 
freedom relevant to real-world inflation, in which case the specific predictions we derive could 
help describe reality. 

Recently, the idea that the Standard Model Higgs field, non-minimally coupled to gravity, can 
lead to infiation was proposed in Ref. [17j. Those authors argued that the radiative corrections 
to the potential are negligible and hence the inflationary parameters can be computed using the 
classical Lagrangian. They found that the cosmological predictions are in good agreement with 
cosmological data, independent of the Standard Model parameters, such as A. On the other 
hand the authors of Ref. |18j criticized their approach, suggesting that the quantum corrections 
to the potential can be very important. They concluded that a Higgs lighter than 230 GeV 
cannot serve as the inflaton, because the predicted spectral index is ruled out by WMAP5 data 
[19j . Ref. |18j only incorporated quantum corrections at leading log order, extrapolated from low 
energies. Here, in contrast, we will compute the full renormalization group improved effective 
action at 2-loops. We conclude that running inflation based upon a Standard Model Higgs makes 
predictions that are consistent with current cosmological data, and leads to firm predictions for 
the PLANCK satellite and the LHC. Our main result is a correlation between the spectral index 
and the Higgs mass, see Fig. [T] This correlation is absent in the classical theory. The origin 
of the correlation lies in the interactions of the Standard Model, which dictate the form of the 
effective action. 

In Section |2] we review inflation with non-minimally coupled scalars. In Sectionjsjwe investi- 
gate the classical theory of the Higgs non-minimally coupled to gravity. In Section mwe describe 
our method for obtaining the quantum corrected effective action. We compute all the inflation- 
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Figure 1: The spectral index rig as a function of the Higgs mass ruf^ for a range of light Higgs masses. The 3 
curves correspond to 3 different values of the top mass: mt = 169 GeV (red curve), mt = 171 GeV (blue curve), 
and mt = 173 GeV (orange curve). The solid curves are for as{mz) = 0.1176, while for mt = 171 GeV (blue 
curve) we have also indicated the 2-sigma spread in asijnz) = 0.1176 it 0.0020, where the dotted (dot-dashed) 
curve corresponds to smaller (larger) as- The horizontal dashed green curve, with ris ~ 0.968, is the classical 
result. The yellow rectangle indicates the expected accuracy of PLANCK in measuring Us (Aris ~ 0.004) and the 
LHC in measuring rrt^ (Amij ft: 0.2 GeV). In this plot wc have set A^e = 60. 



ary observables numerically and present results in Section |5] Finally, we review our results and 
discuss their significance in Section [6j 

2. Non-Minimal Inflation 

Here we briefly review the recipe to compute inflationary observables, which will be used in 
the later sections, and the latest observational constraints. 

Consider a real scalar field non-minimally coupled to gravity via the Ricci scalar TZ. The 
class of effective actions we consider is 



d x\J —g 



\^raij{m-\mm''-y{<i>) 



(2) 



where we allow for a general coefficient of the Ricci scalar /(</>), general coefficient of kinetic energy 
and general potential V{(^\ Here mpi ~ 2.43 x 10^^ GeV is the reduced Planck mass; we 
are effectively assuming that the field is stabilized at the end of inflation with /((/>o) ~ 1, as 
will be the case for the Standard Model Higgs. 

The cosmology of this theory is most easily studied by performing a conformal transformation 
to the so-called "Einstein frame" where the gravity sector is canonical ^rripJZE- This is achieved 



by defining the Einstein metric as gjf^, = f{(t>)gf_iu- The corresponding Einstein frame potential is 
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Furthermore, the kinetic energy in the Einstein frame can be made canonical with respect to a 
new field tr, defined through the equation 



da\ _ k{(l)) 



f{<t>) 



3 2/W 



(4) 



(the second term here comes from transforming the Ricci scalar), hi this frame, the action takes 
the canonical form 



(5) 



which is amenable to straightforward analysis. 

The inflationary dynamics and cosmological predictions is determined by the shape of the 
potential Ve- In the usual way, we introduce the first and second slow-roll parameters, which 
control the first and second derivatives of the potential, respectively. Using the chain rule, these 
are 
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(7) 



where a prime denotes a derivative with respect to 0. Similarly, the third slow-roll parameter C 
is related to the third derivative of the potential as ~ ■m'^^{d^VE / da^){dVE / da) /V'^. 
The number of e-foldings of slow-roll infiation is given by an integral over cj): 
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(8) 



where 



is the value of the field at the end of infiation, defined by e ~ 1. The number of 



e-foldings must be matched to the appropriate normalization of the data set and the cosmic 
history, with a typical value being N^. ~ 60; we return to this point in Section |5] 

The amplitude of density perturbations in fc-space is specified by the power spectrum: 



Ps{k) = Al 



where is the amplitude at some "pivot point" k* , predicted by infiation to be 



Ve 



(9) 



(10) 



and measured by WMAP5 to be A'^ = (2.445 ± 0.096) x lO'^ at k* = 0.002 Mpc"^ [TO]. The 
corresponding spectral index Us, running of the spectral index a = dus/dlnk, and tensor to 
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scalar ratio r, are given to good approximation by 



Us = l-6e + 277, (11) 
a = -246^ + 16e?7 - 2C^ (12) 
r = 16e. (13) 

The combined WMAP5 plus baryon- acoustic-oscillations (BAO) and supernovae (SN) data con- 
siderably constrain and r. Assuming negligible a, as will be the case for running inflation, 
the constraints are: 0.93 < < 0.99 and r < 0.22 (at 95% confidence level). 



3. Classical Analysis 

Without essential loss we can rotate the Higgs doublet so that it takes the form ~ 
(1/V2)(0,u + (p). Only the real field </> will play a role in our analysis. Specializing to gauge 
invariant, dimension < 4 operators, without higher derivatives, the functions /((/)), k((j)), and 
V{(p) must take the form 

/(0)-l + ^, M0) = 1, V{cl>)^^{^'^v')\ (14) 
m^j 4 

where v ~ 246.2 GeV is the vacuum expectation value for the Higgs field, setting the electroweak 
scale. The self coupling A is in one-to-one correspondence with the Higgs mass, namely = 
2Av^. Current experimental bounds on the Higgs mass (and hence A) are as follows: 

1 14.4 GeV < TTih < 182 GeV, 

O.IK A <0.27, (15) 

where the lower bound comes from direct searches and the upper bound comes from a global fit 
to precision electroweak data (95% CL) |20j. 

In this theory, inflation takes place at energies many orders of magnitude above the elec- 
troweak scale ((/)^ ^> v'^). Hence, during inflation the potential is well approximated by the 
quartic potential: V{(f>) = and this form of the classical potential will be sufficient through- 
out this Letter. The corresponding potential in the Einstein frame is then 

VEi^)= (16) 

(1 + ^)^ 

which approaches a constant Vq = ATOp,/4^^ at large field values (p 3> Wpi/v^ (we assume 
^ > 0). This fact allows slow-roll infiation to take place [121 HH HZ]- It is notable that through 
this mechanism slow-roll inflation emerges unusually "naturally" . 

It is useful to define the dimensionless quantity ip = -^f 0/mpi which controls the cosmological 
evolution: infiationary stage {ip 3> 1), the end of infiation (-0 ^ 1), and the low-energy regime 
{tp ^ 1). Indeed the potential Ve plotted in Fig. |2] displays the familiar quartic behavior for 
small ip values, but asymptotes to a constant for large ^. 
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Figure 2: The potential in the Einstein frame Ve, normalized to a reference value Vq = Amp[/4g^, as a function 
of the Higgs field 1/1 = \^ <j>/mpi. The dashed green curve is the classical case (independent of Higgs mass), 
the solid blue (red) curve is the quantum case with Higgs mass mf^ = 126.5 GeV (mji = 128 GeV). We have set 
rrit = 171 GeV and as{mz) = 0.1176 for this plot. The inset focusses on the slow-roll inflationary regime. 



Using eqs. Q, ([6]), and the slow-roll parameters are readily computed. The exact results 
are not very transparent. They simplify for large ^, which is the case of physical interest: 
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(17) 



We see that at large ip (during slow-roll inflation) 77 is dominant, and will primarily control the 
predictions for the spectral index. The number of e-foldings is computed from eq. ([8|) giving 



1 + V^^ 

1 + 



2 

end 



(18) 



where ^Acnd — (4/3)^^'' is the value of ip at the end of inflation (e ~ 1). Eqs. ( 17 1 and ( 18 1 provide 
a parametric description of e(A^e), vi^e), and C(^e), thus determining Ug, a, and r as a function 
of Ne, i.e., we can trade the unknown value of the Higgs fleld during inflation (f){— rripiip / \/^) for 
the number of e-foldings N^. 

For Ne — 60 we find the following results for the spectral index, the running of the spectral 
index, and the tensor to scalar ratio: 



0.968, 



-5.2 X 10" 



3.0 X 10-^ 



(19) 



We see that a and r are rather small. This will remain qualitatively true in the quantum theory, 
but the corrections to ns are quite important, as we explore in detail in the next section. 

1/^ 



/Ne, the amplitude of 



Finally, using eq. ( 10 1 and expanding to leading order in 1 /ip 
density fluctuations is found to be 



^2 727r2 
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Since this must be 0(10"^), it is impossible to satisfy for A = 0(0.1) and ^ = 0(0.1) (which 
might be considered "natural" values). One possibility is that A is extremely small, but that is 
incompatible with experimental bounds on the Higgs mass, see eq. ( [Tsl ), and is not stable under 
renormalization. Instead, following [17', we assume ^ = 0(10*) in order to obtain the correct 
amplitude of density fluctuations with A = 0(0. 1) . The need to dial a parameter to large or small 
values, so that is consistent with observations, is a common feature to all known inflation 
models. It will also apply in the quantum theory. 



4. Quantum Analysis 

We now consider how quantum corrections modify the classical results of the previous section. 
In order to do so, we need to compute the effective action that takes into account the effects of 
particles of the Standard Model interacting with the Higgs boson through quantum loops. The 
frame we calculate in is the original "Jordan" frame which defines the theory. The quantum 



theory modifies all three functions /(</>), k{<j)), y(0) from the classical expressions in eq. (14 1. 

The quantum corrections to the classical kinetic sector k((f)) — 1 arise from wave-function 
renormalization, and are approximately ^-independent. It is simple to check that at large ^ the 
second term in eq. Q scales as ^ 1, while the contribution from the k{(j))/f{4>) term scales 
as 1/^. Hence corrections to k((j)) occur with a factor 1/^, in addition to suppression by loop 
factors and couplings. 

The quantum corrections to the classical gravity sector /(0) = 1 + are more subtle. 

Let us start by considering the case of a (classical) background gravitational field. In this case 
the conformal anomaly induces a 1-loop /3-function for ^ given by (21j 



2X + y^-lg^~y^ 



(21) 



The term proportional to A, coming from Higgs running in a loop (see Fig. [sj^a)), is potentially 
important during inflation. We will return to this point soon when we include the (classical) 
back reaction of gravity, and argue that in fact this contribution is negligible. The remaining 
terms arise from external leg corrections and cancel against wave-function renormalization to 
good approximation. Hence corrections to /((/>) are ignorable also. 

Finally we turn to the computation of the potential sector V{(p). Let us begin with the 
flat space analysis. The RG improved potential for the Higgs in the Standard Model is (see 
e.g. Ref. [22] for a review) 

v{c^)^\\m)Gmr^\ (22) 

{(j) ^^) where t{(j)) = ln((/)//x), and fj, is the normalization point; taken to he ii — mt in this 
Letter. Here X{t) encodes the running of A, while G{t) = exp(— dt'^{t')/{l + j{t'))), where 
7 is the anomalous dimension of the Higgs fleld, encodes wave-function renormalization. The 
running of A is governed by the renormalization group equation: dX/dt = Px/{1 + 7). At 1-loop 
it is 



/3a 



(4^)2 



24A^ -6yt+ + -/g'^ + + KUyj " 9.g^ - 3/^) 



(23) 
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Figure 3: Some representative Feynman diagrams. Top row: renormalization of the conformal coupling ^ with 
Higgs in loop (a), and renormalization of top quark's Yukawa coupling with gauge boson (b) and Higgs (c) across 
vertex. Bottom row: renormalization of quartic coupling A with Higgs (d), top quark (e), and gauge boson (f) in 
loop. 



At low energies, the two most important terms here are the self coupling 24 (see Fig. [sj^d)), 
which tries to drive A to large positive values, and the top quark —6 yf (see Fig. [sje)), which tries 
to drive A towards zero. This is summarized in Fig. |4j This leads to a delicate interplay between 
the Higgs mass and the top mass. For irih ^ mt, the 24 A^ term dominates and A will eventually 
hit a Landau pole at high energies. For nih ^ nit, the —Qyf dominates and A will go negative 
which is a sign of vacuum instability. The "Goldilocks" window for the Higgs mass, where the 
theory is both perturbative and stable up to very high energies is also the regime in which the 
quantum corrections are relatively small, allowing for slow-roll inflation. At high energies, the 
contribution from gauge bosons (see Fig.jsjf)) are important and increase A. 

In the recent work of Barvinsky et al. |18j the top quark's Yukawa coupling was approximated 
by the tree level value: yt = y/2mt/v for all energy scales. This provides a significant negative 
contribution to (3\, forcing A to negative values and vacuum instability in large regions of pa- 
rameter space. Instead it is essential to include the running of the top Yukawa coupling in the 
analysis: 

'^^^ " (4^) 

which is negative due to the large negative contribution from the strong coupling —Sg^yt (see 
Fig. |3jb)). Hence yt runs to smaller values at high energies; see Fig. |4] 

In our work, we have included the complete running of the 5 couplings: A, yt, gs, g, and 
g' to 2 loops, to ensure accurate resultsj^ The /3-functions are summarized in Appendix |a] 
Furthermore, we have adopted the pole mass matching scheme for the Higgs and top masses. 



9 2 o 2 9 2 1'' /2 



(24) 



"The 3-loop running is unknown for the Standard Model, but would need to be abnormally large to have an 
effect. 
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Figure 4: This plot summarizes some of the most important effects of the renormahzation group fiow. The 
red curve shows the running of the quartic coupUng A(t)/A(0) for a light Higgs rrif^ = 126.5 GeV. The dotted 
purple curve is the top running yt{t)/yt{0) Eind the dot-dashed cyan curve is the commutator function s{t), with 
^ = 2.3 X 10^ and fi = mt- The right-hand region is the slow-roll inflationary regime; here A rises (and so ris does 
too), as highlighted by the inset. 



given in the Appendix of [23' . For the sake of brevity, we do not reproduce the pole matching 
details here. 

We now consider the effective potential V{(f)) including the effect of the non-minimal cou- 
pling to gravity £^(j)^TZ. The calculation is difficult to perform exactly. However, we can obtain 
approximate results for large ^ fairly simply. Following [12 , one can heuristically identify a non- 
standard commutator for (j) as follows. From eqs. Q and Q we see that when the gravity sector 



is canonical, the kinetic sector is non-canonical 
the canonical momentum corresponding to (j) is 



-\{dE4>Y {du/dcj)) . On a spatial hypersurface. 



dL 



V-9E {9[ 



da 



da\- 
d^J 



(25) 



where is a unit timelike vector. Imposing standard commutation relations for (f> and tt, we 
learn that [^(x),0(y)] = i h s{(j)) 5^^\:x. - y) , with 



da 



1 



1 + ^ 
(6^ + 1)5^ 



(26) 



For (f) ?Tipi/^ (the low energy regime) we recover the ordinary value of the commutator s = 1, 
while for (j> ^ '7ipi/\/^ (the inflationary regime) we see a suppression in the commutator by a 
factor of s = 1/(6^ + !)• So in the inflationary regime with ^ ^ 1, quantum loops involving the 
Higgs field are heavily suppressed. 

To summarize, our prescription for the renormahzation group improved effective potential in 
the presence of non-minimal coupling is to assign one factor of s{(f>) = s{fie^) for every off-shell 
Higgs that runs in a quantum loop. This factor is plotted as the dot-dashed cyan curve in Fig. |4] 
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In eq. (23 1, for example, this prescription means the replacement 24 A 24: s A , as that term 
arises from two Higgs ofF-shell propagators, while all other terms are untouched since they only 
involve other fields in loops (see Appendix |A] for more details) . This provides an important 
modification to the high energy running of couplings, and explains why the running of ^ from 
the diagram of Fig. [sj^a) is suppressed. Apart from this modification, the RG improved analysis 



is as standard, as summarized in eq. (22 1. We have checked our prescription against detailed 
analytical calculations of the effective action of non-minimally coupled scalars in the literature 
(e.g., see [24l |25]) and have found excellent agreement. We assume that quantum corrections 
from graviton exchange are small, see Appendix |B] 

5. Results and Predictions 

After numerically solving the set of 5 coupled renormalization group differential equations of 
Appendix [A| for the couplings: A, yt, gs, <?, and g' , we have obtained the effective potential V(0) 
in the full quantum theory, as a function of input parameters, such as the Higgs mass. Some 
representative potentials in the Einstein frame are given in Fig. (|2|. The inset clearly exhibits 
variation of the effective potential Ve{(I)) with Higgs mass, which was absent in classical case. 
As we lower the Higgs mass, approaching the instability, the magnitude of the first derivative 
is raised and the that of the second is lowered (see the blue and red curves). This leads to 
modifications to the cosmological parameters. 

Following the recipe we outlined earlier in Section [2] we are able to efficiently compute the 
spectral index in the RG improved theory using Mathematica. Recall that in the classical theory. 
Us is independent of the parameters of the Standard Model, and its value was found to be 
Us — 0.968 (for = 60). In the quantum theory, we find that Ug depends on several of the 
Standard Model parameters, in particular on the Higgs and top masses, see Fig. [T] As the top 
mass is varied through its experimentally allowed range (169 GeV ^ nit < 173 GeV) the spectral 
index varies noticeably. In particular, as we lower the Higgs mass towards vacuum instability, the 
spectral index increases substantially. To achieve successful infiation with < 0.99, we require 

„.,. > .„G.V.3.G.V (^5^^) - ( "'too'.:"" ) (-) 

where (5 ~ 2 GeV indicates theoretical uncertainty from higher order corrections (such as 3-loop) . 
This bound almost coincides with that from absolute stability presented in Ref. [23 . Note that 
near the boundary A is small, so the corresponding ^ to obtain the observed is reduced from 
its classical value ^ ~ 10'' by an order of magnitude or so to ^ ^ 10'^. 

Let us now trace the chain of logic behind the rise in Us- For a light Higgs, P\ is dominated by 
the top and gauge boson contributions. For a heavy top, the top contribution is dominant at low 
energies, causing f3\ to be negative and thus driving A to low values as the energy is increased. 
At the same time, the top Yukawa coupling runs, with dominant contributions coming from 
gauge fields and Higgs running in a loop, with the gauge fields slightly dominant causing yt to 
decrease with energyj^ At very high energies (p ^ nipi/y/^ (the infiationary regime), the Higgs 
running in the loop is highly suppressed, causing yt to jump to even lower values. Hence the 



^Note that the 2-loop term — lOSytg^ /(47r)'' in I3y^ (see eq. (kMI) speeds up the running compared to 1-loop. 
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top contribution to the running of A becomes subdominant, the gauge boson contributions now 
dominate and A rises, as seen in Fig. |4] (inset). Since A is concave up, this increases 77 and hence 
the spectral index. 

In Fig. [1] and in all plots we have chosen the reference value A^o = 60. For N,, close to Nq, we 
can Taylor expand rig to linear order: 

(ITl 

ns{N,) = n,(7Vo) + ^(iVe - iVo) + • ■ • ■ (28) 

Now, the spectral index is in fact a function of all the parameters, including and ^: Us — 
ns(iVe, ^, . . .). As in the classical theory, we have fixed ^ such that the amplitude of density 
fluctuations is in agreement with observations (requiring ^ ^ 10^\/A). In this way, we can think 
of ^ = ^(A^'e), so from the chain rule 

dris dus dus rf^ 



(29) 



The first term is precisely the (negative) of the running of the spectral index a = dris/dlnk, 
while the second term is found to be very small numerically. Hence to a good approximation we 
can write 

n,{Ne)~ns{No)-aiNo){N,-No). (30) 

We plot the running of the spectral index a{No — 60) in Fig. [s] (left). We see that a « —5 x 10^*^ 
(as in the classical case), with some variation for low Higgs masses as we approach the instability. 
However, this is still far too small to be detected by PLANCK, which is expected to be only 
sensitive to a = 0(10"^) [5^. Hence the main usefulness of Fig. [s] (left) is that it should be used 
in accompaniment with Fig. [T] and eq. ( 30 1 to infer the value of Us for different values Ne (as 



long as Ne does not vary too far from A^o = 60). 

The actual number of e-foldings of inflation is related to the wavenumber of interest k, the 
energy density during inflation Ve, the energy density at the end of inflation Knd, and the energy 
density at the end of reheating [B] 

jVe-62-ln— ^ + -ln , + - In — - — In — . (31) 

aoHo 4 (10i6GeV)4 4 K„d 12 ^ ' 

Since the Higgs is strongly coupled to Standard Model fields, reheating is expected to occur 
automatically. As has only a weak dependence on p^Gh, the details of reheating are rather 
inconsequential to our mass bounds, but may be calculable ^7\. According to [28], r„h ~ 
-[^q13.5 QgY giving — 59 for the classical theory. In our case, we must take into account the 
variation in the scale of inflation due to the quantum corrections. In Fig. [S] (right) we plot r 
versus the Higgs mass. Since we have fixed ^ such that the amplitude of density fluctuations 
is at the observed value, the energy density of inflation Ve is simply proportional to r. Using 



eqs. (10 1 and (13), we have 

Ve = 2-'< r « (7.8 x 10^^ GeV)^ (^) . (32) 
Since r changes by a factor of order 2, as we vary the Higgs mass, then Ve changes by the same 



amount. From eq. (31 1, rescaling Ve and V^^^ by a factor of p, say, the number of e-foldings is 
shifted by AiVg = ^ logp, which is « 0.1 for p = 2. Hence the variation in with the Higgs 
mass is very small. 
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Figure 5: The running of the spectral index a(xlO*) (left panel) and the tensor to scalar ratio r(xlO^) (right 
panel) as a function of the Higgs mass m,fi- The 3 solid curves correspond to 3 different values of the top mass: 
mt = 169 GeV (red curve), mt = 171 GeV (blue curve), and mt = 173 GeV (orange curve). The horizontal dashed 
green curve, with a ~ —5.2 X 10"* and r ~ 3.0 X 10~^, is the classical result. We have set = 0.1176 and 
A^e = 60 in this plot. 



6. Discussion 

A number of papers have discussed bounds on the Higgs mass coming from demanding stabil- 
ity of the vacuum, e.g., see [231 [251 EH] ■ Cosmological constraints only require metastability on 
the lifetime of the universe, which places the constraint ruh > 105 GeV f23|. However, if we further 
demand that the Higgs drive inflation, we find that heavier Higgs are required: ruh ^ 126 GeV 



(depending on the top mass, see eq. (27 1), which essentially coincides with the bounds from abso- 
lute stability. Furthermore, by demanding that the theory remains perturbative to high energies 
{nih ^ 190 GeV), we establish a correlation between both stability and triviality bounds, and 
inflation. 

More precisely, we have established a mapping between the renormalization group flow and 
the cosmological spectral index. Over a substantial range of parameter space the classical value 
rig ~ 0.968 (for N^, — 60) emerges as a good approximation, but there are corrections. Given 
a detailed microphysical theory, such as the Standard Model, we can explicitly calculate such 
corrections, as summarized in Fig. [T] This plot displays a sharp rise in the spectral index towards 
0.98, or so, as we approach vacuum instability for a light Higgs. 

It is likely that the Standard Model is only the low-energy limit of a more complete theory, 
accommodating the facts that do not flnd explanation within the Standard Model, such as neu- 
trino masses, dark matter, baryon asymmetry, etc. Our methodology is still applicable, so long 
as we can control the relevant /3-functions. 

In principle some quite different scalar fleld, not connected to the Standard Model Higgs, 
could drive running inflation. The central requirement is a large coefficient for the 0^7?. term. 
It is possible that such a coefficient could emerge as some sort of Clebsch-Gordan coefficient, or 
from the coherent addition of several smaller terms (involving more basic scalars (fij). It is also 
possible to consider, in the same spirit, the dimension 3 interaction (j/IZ, which arises for generic 
scalar fields, though not of course for the Standard Model Higgs. Furthermore, as discussed in 
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Appendix [B] and Refs. [351 ISS]) the inclusion of higher dimension operators may significantly 
affect the predictions of the original theory and even spoil its validity (as is the case in many 
inflationary models). These possibilities, and their possible embedding in unified field theory or 
string theory, deserve further investigation. 

If the Higgs boson exists and it is in the mass range considered in this Letter, the LHC will 
discover it and will determine the Higgs mass with a precision of about 0.1% [21], which means an 
uncertainty Am/j « 0.2 GeV. In order to extract accurate correlations between the inflationary 
observables and the Higgs mass it is crucial to improve the precision with which we know the other 
parameters of the Standard Model, in particular the top quark mass and the strong coupling. 
The current value of the top mass from direct observation of events is = (171.2 ± 2.1) GeV 
[5D]. In the near future, the LHC will improve the determination of the top mass, but relatively 
large systematic uncertainties will prevent a top mass determination to better than 1 GeV; more 
conservatively, the top mass will be determined at LHC with an error Amt = l-;-2 GeV. Looking 
further ahead, the ILC is expected to be able to measure the top mass to ~ 100 MeV. So together 
with the measured Higgs mass from the LHC and improved precision on the strong coupling, as 
well as calculating higher order effects and reheating details, running inflation in the Standard 
Model will predict a rather precise value for the spectral index. 

We would like to thank M. Amin, F. D' Eramo, A. Guth, and C. Santana for useful discussions, 
and A. Riotto for comments on the manuscript. We thank F. L. Bezrukov, A. Magnin, and 
M. Shaposhnikov for correspondence. The work of ADS is supported in part by the INFN 
"Bruno Rossi" Fellowship. The work of ADS, MPH, and FW is supported in part by the U.S. 
Department of Energy (DoE) under contract No. DE-FG02-05ER41360. 

Note Added 

Our Letter appeared simultaneously on the arXiv with Ref. [32 , which also studied the 
quantum corrections to inflation driven by the Standard Model Higgs. The central conclusion 
of both papers is that the classical analysis provides a good approximation over a wide range of 
parameters, but that quantum corrections are calculable and can be quantitatively signiflcant. 

For a top mass of nit = 171.2 GeV, Ref. ^32] found that in order to have successful inflation 
the Higgs mass is constrained to be in the range: 136.7 GeV < m/j < 184.5 GeV, and the spectral 
index decreases from its classical value as rufi approaches the lower boundary. In this note we 
briefly discuss the similarities and differences between their analysis and oursj^ 

In our analysis, we computed the full RG improved effective potential. We did this including 
(i) 2-loop beta functions, (ii) the effect of curvature in the RG equations (through the function 
s), (iii) wave-function rcnormalization, and (iv) accurate specification of the initial conditions 
through proper pole matching. On the other hand, f32| did not compute the full effective potential 
or include any of the items (i)-(iv)j^ Instead Ref. 32J approximated the potential at leading log 
order with couplings evaluated at an inflationary scale after running them at 1-loop (this is one 
step beyond [18j where couplings were not run). 



*The discussion here refers to version 1 of I32| . 

^Though wave-function renormalization was not included in |32l , external leg corrections in the running of A 
were included. However these two effects roughly cancel against one another. 
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The lower bound on the Higgs mass we find in eq. (27) is about llGeV lower than that 
found in Ref. (m/j > 136.7 GeV). This numerical discrepancy is due to several of the above 
simplifications, but the dominant difference comes from (i) inclusion or not of 2-loop effects 
(importantly, the — lOSytgf /(47r)* term in fiy^, see eq. (34l), and a second significant difference 
comes from (iv) pole matching. Higher order effects (such as 3-loop) and uncertainty in the 
strong coupling as{mz) also modify the bound, as we summarized in eq. (27 1. 

A precise upper bound on the Higgs mass {nih < 184.5 GeV) is stated in [32]. The basis of 
this is the famous "triviality bound", see e.g., Ref. [33j, which has little to do with inflation. 
The theory ultimately requires a cutoff, and exactly how low a cutoff one feels comfortable with 
(or equivalently, how large a value of A one regards as acceptable) is arguable. We feel that our 
stated semi-quantitative bound nih ^ 190 GeV adequately represents the situation. 

In |32l ris decreases as m/j approaches its minimum value, while we find that increases 
(see Fig. [T|). This behavior depends critically on the value of yt during infiation, as compared 
to the value of the gauge couplings. If yt is small, then Ug increases, and vice versa. Ref. [32] 
overestimated yt during inflation and hence obtained the opposite behavior. This is primarily 
due to ignoring items (i) and (ii) above. Ignoring (i) misses the 2-loop term — 108^4(7^ /(47r)'' in 
/3yj, and ignoring (ii) maintains the 1-loop term |yj'/(47r)2 in f3y^ during inflation. 

Finally, |32j computes quantum corrections with both a field-independent cutoff (as we use) 
and a field-dependent cutoff in the original "Jordan" frame. Either procedure defines a possible 
model, but the field-independent cutoff is more in the spirit of the motivating arguments, based 
on dimension < 4 effective Lagrangians. 



A. 2-Loop RG Equations 



In this appendix we list the RG equations for the couplings A, yt, g' , g, gs at energies above nit 
at 2-loop In each case, we write dX/dt — f3\/{l + 7], etc., where t ~ ln(j)/fi. Also, we insert 
one factor of the commutator function s{iJ,e*) (see eq. (26 1) for each off-shell Higgs propagator p°] 

For the Higgs quartic coupling we have 



/3a 



1 



(47r)2 
1 

(47r)4 
A 



1 

48 



6yt + I (25^ + (, 
(915/ - 289/g'2 „ 559g^g'^ _ 379g'^) + SOsy 



9^2 _ 3^.2 ^ ^2y2) A 
Sg" , 



Vt 



i2gl + 3sA 



73 ^ 39 2 



629 



Vt 



4 
21 



g g 



n , --„„,4 , io8sVA + 365V'A-312s4a2 



24 



sg 



^.'^ + A (f,2 + f,'2 + 80,2 -144.2a)) 



(33) 



We have carefully extracted out all Higgs propagators contributions at 1-loop order by the appropriate insertion 
of factors of s. For the 2-loop contributions we have only inserted s for the obvious terms. The complete set of 
insertions are tedious and provide negligible corrections. 
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For the top Yukawa coupling we have 



yt 



■'Vt 



(47r)2 
19 



Vt 



(47r)4 



225 

-16' 



131 



For the gauge couphngs gi ~ {9',g,9s} we have 



1 

(4^) 



3 

4 



(34) 



(35) 



ith 



199/18 9/2 44/3 
((40 + s)/6,-(20-s)/6,-7), B=\ 3/2 35/6 12 



rf* = (17/6,3/2,2). (36) 



11/6 9/2 

Finally, the anomalous dimension of the Higgs field is 



-26 



7 



1 

" (47r)2 
1 

(47r)4 



^ 3.9'2 „ 2 



271 



9 
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5 /9 



2 V4 



17 



9' + j^9'' + Sg',]yt+^^yt-Gs'^ 



3 \2 



(37) 



B. Remarks on Running Inflation as an EFT 

The Lagrangian analyzed in this Letter is not renormalizable in the conventional sense, nor 
is it "technically natural" from the point of view of effective field theory. In this appendix we 
remark on the validity of such a theory at high energies (for related discussions see [35l|36]) and 
elaborate on the spirit of our calculations. 

The novelty of running inflation is to introduce the non-minimal coupling ^0^7?. into the low 
energy Lagrangian, which is allowed by all known symmetries of the Standard Model and gravity. 
This term is dimension 4 in the same sense that the kinetic term g'^'^ d^(j)du4> is also. However, if 
we expand around flat space g^^, = rj^^ + /i/iiy/n^pi then the new term is dimension 5 at leading 
order, plus an infinite tower of corrections 

^(j)'^n - ^(/i^n/i/mp, + 

which is connected to the non-renormalizability of gravity in 4 dimensions. This suggests that 
non-minimally coupled theories becomes strongly interacting at scales A ^ mpi/^. This can be 
compared to minimally coupled theories with A ~ mpi. 

Without any protecting symmetry, we cannot forbid infinite towers of corrections to the 
dimension 4 effective Lagrangian £4, including those of the form 

ri>0 ^ ^ n>0 
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which applies to non-minimal models (with A ~ rripi/^) and minimal models (with A ^ m-pi). The 
values of the higher order Wilson coefficients a„, 6„ cannot be determined without knowledge of 
the behavior of gravity at energy scales above A, since these terms arise from graviton exchange. 
If we take a naive estimate o„, 6„ = 0(1), then the required flatness of the inflationary potential 
is jeopardized. This applies both to running inflation and to many minimal inflation models, such 
as m^(^^ chaotic inflation, since in both cases cj) > A during inflation. As there is no increased 
symmetry in the limit an,bn 0, such theories are not "technically natural". 

On the other hand, we currently have no evidence for a^, = C(l), as these terms arise from 
graviton exchange, whose effects are yet to be seen in any experiment. There does exist a logical 
possibility that graviton exchange at high scales is softer than naive estimates suggest (Ref. [37] 
may be an example), rendering a„,&„ small, preserving unitarity, and leaving our calculated 
potential V{(j)) essentially unaltered. It is in this spirit of including only the known Standard 
Model loops, and not those of unknown graviton loops, that we have obtained our results and 
predictions - which are highly falsifiablc. 
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